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Oh" Abstract 

T S^ , We propose an A4 flavor-symmetric SU(3) ® SU(2) <g> U(l) seesaw model where lepton num- 

ber is broken spontaneously. A consistent two-zero texture pattern of neutrino masses and 

(N ■ 

£> ■ mixing emerges from the interplay of type-I and type-II seesaw contributions, with important 

00 ' 

q^ I phenomenological predictions. We show that, if the Majoron becomes massive, such seesaw sce- 

00 ' 

o 



nario provides a viable candidate for decaying dark matter, consistent with cosmic microwave 



[~^. . background lifetime constraints that follow from current WMAP observations. We also calculate 

the sub-leading one-loop-induced decay into photons which leads to a mono-energetic emission 



line that may be observed in future X-ray missions such as Xenia. 
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INTRODUCTION 



The discovery of neutrino oscillations [1H5J, now confirmed at reactors and accelera- 



tors [6|-|8|, has brought neutrino physics to the center of particle physics research. Global 
analysis of current oscillation data indicate that the pattern of lepton mixing differs 
sharply from that characterizing quarks [9[. Understanding the origin of neutrino mass 
and the pattern of neutrino mixing angles from basic principles constitutes a major chal- 
lenge 10|, |ll|. A paradigm framework to generate neutrino masses is p rovided by the 
seesaw mechanism, for which several realizations have been proposed [12[. The observed 
pattern of neutrino mixing may arise from suitable non-abelian flavour symmetries, as 



those based on the A4 group |13l— Il6| . 



Elucidating the nature of dark matter constitutes another intriguing problem of mod- 
ern physics which has so far defied all efforts. It is therefore crucial to build a fundamental 
particle physics theory of dark matter and, since the Standard Model of elementary parti- 
cles (SM) fails to provide a dark matter candidate, such theory necessarily requires physics 
beyond the SM. 



Here we suggest a version of the seesaw mechanism containing both type-I 17H24 1 



and type-II contributions 



23 



in which we implement an A4 flavor symmetry with 



22 



24J | . We study the resulting pattern of vacuum 



spontaneous violation of lepton number 

expectation values (vevs) and show that the model reproduces the phenomenologically 

consistent and predictive two- zero texture proposed in Ref. [29 1. 

In the presence of explicit global symmetry breaking effects, as might follow from grav- 
itational interactions, the resulting pseudo-Goldstone boson - Majoron - may constitute 
a viable candidate for decaying dark matter if it acquires mass in the keV-MeV range. 
Indeed, this is not in conflict with the lifetime constraints which follow from current cos- 
mic microwave background (CMB) observations provided by the Wilkinson Microwave 
Anisotropy Probe (WMAP) 30[. We also show how the corresponding mono-energetic 
emission line arising from the sub-leading one-loop induced electromagnetic decay of the 
Majoron may be observed in future X-ray missions [311 ] . 

The paper is organized as follows. In section [XT] we describe our A4 model while in 
section IIIII we discuss the symmetry breaking structure which is required to obtain the 
correct neutrino texture. In section HVl we update the neutrino parameter analysis and we 
study the implications of a decaying Majoron dark matter scenario in section [V] Further 
discussion is presented in the concluding section IVII 



II. THE MODEL 

Our model is described by the multiplet content specified in Table |T] where the trans- 
formation properties under the SM and A 4 groups are shown (as well as the corresponding 
lepton number L). The Li and lm fields are the usual SM lepton doublets and singlets and 
v R the right-handed neutrinos. The scalar sector contains an SU(2) triplet A, three Higgs 
doublets $,/ (which transform as a triplet of At) and a scalar singlet a . Three additional 
fermion singlets Si are also included. 
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Taking into account the information displayed in Table HJ and imposing lepton number 
conservation, the Lagrangian responsible for neutrino masses reads 

-C L = hxLx {y R $)' x + h 2 L 2 {v R $>) 1 + h 3 L 3 {y R $){ 
+ \L{CAL 2 + \L T 2 CAL X + X'L^CAL 3 
+ M R {S~ L y R ) l + h (SlCS^ a + h.c. , (1) 

where h and A are adimensional couplings, M R is a mass scale and 

^ -A+/V2 A« )'*' [ft)- [I> 

Note that the term (u R Ci' R y i a is allowed by the imposed symmetry. This term however 
does not contribute to the light neutrino masses to the leading order in the seesaw ex- 
pansion and we omit it. Alternatively, such term may be forbidden by holomorphy in a 
supersymmetric framework with the following superpotential terms 

W = • • • + \e ab h»L*i> c H b u + M R u c S + \hSSa 
where the hats denote superfields and the last term replaces the corresponding bilinear 



employed in Ref. 32|, l33(. Assuming that the Higgs bosons $j, A and a acquire the 



following vevs (see section II III below) 

<0?> = (4) = (4) 



x/3' 



(A ) = « A , {a} = u a 



(3) 



we obtain an extended seesaw neutrino mass matrix M 32|-|34J in the (ul, v c , S) basis 



/ 
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m D = v diagOi, /i 2 , h) U, U = — = 



1 1 1 

1 u; u 2 



(4) 



with a; = e 2m ^ 3 , M = Mr diag(l, 1, 1) and \x = u a h diag(l, w 2 , w). This leads to an 
effective light neutrino mass matrix Mr v given by 



A/rT- 1 ., nr-l^T hv U 



Ml = moM 1 fiM-'mj) 



Ml 



( hi \ 

h 2 h 3 
h 2 h 3 



(5) 



On the other hand the vev of the triplet, «a, will induce an effective mass matrix for the 
light neutrinos from type-II seesaw mechanism 



/o A 0\ 



M 1 * = 2u A 



A 
A' 



(6) 



and the total effective light neutrino mass matrix will then be 

M v = Ml + M 1 } . 



(7) 



In Ref. 



291 ] it was shown that the neutrino mass matrix given by Eq. (|7j) could explain 



the currently available neutrino data. In section [IV] we will present an update of that 
analysis taking into account the latest neutrino oscillation data. 



III. A A INVARIANT HIGGS POTENTIAL 

We now address the question of the minimization of the neutral Higgs scalar potential, 
which is a necessary condition to reproduce the structure of the neutrino mass matrix 
presented in the previous section. With the assignments of Table [H the Higgs potential 



consistent with gauge and A4 invariance and lepton number conservation reads, 

V = V{$) + V($,A,<t), (8) 

where V(<&) is given as (the decomposition of the tensorial product of two triplets in A4 
is shown in the Appendix): 

V($) = m% ($ t $) 1 + Ax ($ t $) 1 ($ t $) 1 + A 2 ($ t $) 1 , (&$) v , 

+ A 3 (&*) 3s ■ (&*)„ + A 4 ($t$) 3s . ($t $ ) 3Q + As (^t^)^ . ($ t $ ) 3a ; (9) 

and V(<&, A, a) contains pure A , a terms, together with others involving mixed invariant 
combinations of the scalar fields. Assuming the so-called seesaw hierarchy ma C t; < 



u a 24J l , the relevant terms in V(<&, A, a) are 



2 



V($, A, a) = (M 2 A + p\a\ 2 ) Tr(A f A) + A ff |a| 4 + [m 2 a + £ ($ f $) J M 2 - (5$ T A$a* + h.c), 

(10) 
Taking the vacuum alignment for the Higgs doublets $ a given in eq. ([3]) the minimization 
of the Higgs potential with respect to A gives 

AT/ 

_ = =}► (M 2 A + pu^) n A - ^V = . (11) 

We stress that the A4 symmetry, together with the doublet vev alignment assumed in 
Eq. ([3]), requires that the product $tg>$ ~ 1 under v4 4 . If $£g>$ ~ 1', 1", then the second 
term in the above equation would reduce to 25(1 + u + uj 2 )u a = implying u& ~ 0. 
Moreover, as a direct consequence of the requirement $ £g> $ ~ 1 under A4, A and <r must 
have the same (singlet) transformation properties under that group. 
The above equation leads to the following solution for the triplet vev 

5v 2 u a 5v 2 

ua= MiT^r^? (12) 

where the last approximation holds for Ma <C u a . This result shows that the "vev-seesaw" 



1 In contrast to the inverse seesaw models used in Refs. 22b l2J] here we consider large values of u a , 
u a > 10 7 GeV or so. 

2 Notice that the scalar potential contains other invariant terms such as $t$Tr(A^A), Tr(A^A)|er| 2 , 
[Tr(A'l'A)] 2 , etc. Assuming the vev hierarchy ma << v << u a and that the adimcnsional coefficients 
of these terms are of the same order of the ones in V(Q, A, a), then V($,A,er) is enough for our 
purposes. 



relation u&u a ~ v 2 is fulfilled. The minimization with respect to the $ a gives 

W = ° ^ ^$ + 2 ^ M ' ~ 4 ^ m a« ct = 0. (13) 

Finally, 

AT/ 

= o => 2A ff i4 + (m 2 + ^ 2 + pi^) m ct - 2(fo 2 u A = 0. (14) 

0(7 

which, in the limit u&,v << u ff , has the approximate solution 



m? 



"2^' (15) 

as it is typical from spontaneous symmetry breaking scenarios. In summary, we have 
shown that in our framework it is possible to achieve a consistent minimization of the 
scalar potential with non-zero vevs satisfying the "vev-seesaw" relation u/\u a ~ v 2 . 



IV. NEUTRINO PARAMETER ANALYSIS 

Given the two contributions to the light neutrino mass matrix discussed in Eqs. (131) 
and ([6]) one finds that the total neutrino mass matrix has the following structure: 



M v 



(a b 0\ 

6 0c 
\0cdl 



(16) 



This matrix with two- zero texture has been classified as Bl in 35( . One can show that con- 
sidering the (Li, L 2 , L 3 ) transformation properties under A4 as being (1', 1", 1) or (1", 1', 1) 
an effective neutrino mass matrix with .M„(l,2) = 7Vf„(3,3) = is obtained (type B2 



in 35]). Moreover, by choosing A, a ~ V and appropriate transformation properties of 
the Li doublets, we could obtain the textures Bl and B2 as well. Still, the configuration 
A, o ~ 1 would lead to textures which are incompatible with neutrino data since, in this 
case, both type I and type II contributions to the effective neutrino mass matrix would 
have the same form. Since the textures of the type Bl and B2 are very similar in what 
concerns to neutrino parameter predictions, we will restrict our analysis to Bl, shown in 
flU). 

In general, the neutrino mass matrix is described by nine parameters: three masses, 
three mixing angles and three phases (one Dirac + two Majorana). From neutrino oscil- 
lation experiments we have good determinations for two of the mass parameters (mass 
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FIG. 1: Correlation between the neutrinoless double beta decay amplitude parameter \m ee \ and 
the atmospheric mixing parameter. Experimental sensitivities are also given for comparison. 



squared differences) and for two of the mixing angles (0 12 and 6*23) as well as an upper- 
bound on the third mixing angle #13. Using the 3<r allowed ranges for these five parameters 
and the structure of the mass matrix in Eq. (fT6l) we can determine the remaining four 
parameters. The phenomenological implications of this kind of mass matrix have been 
analysed in Refs. [29( and [36|. Here we will update the results in light of the recently 
determined neutrino oscillation parameters [9j. 

The main results are shown in Figs. [1] and [21 In figure [1] we plot the correlation of the 
mass parameter characterizing the neutrinoless double beta decay amplitude: 



m P 



12 2 1 2 2 lia 1 2 2iB I 

\c l3 c 12 m 1 + c 13 s 12 m 2 e + s 13 m 3 e p \ 



(17) 



with the atmospheric mixing angle 6*23. Here Cjj and s^ stand for cos 6ij and sin^ 
respectively. At the zeroth order approximation mi/m 3 = tan 2 #23, and therefore #23 < 45° 
for normal hierarchy (NH), while #23 > 45° for inverted hierarchy (IH). The main result 
from this plot is a lower bound on the effective neutrino mass:|m ee | > 0.03 eV. For 
comparison the range of sensitivities of planned experiments as well as current bounds is 
also given. Note that the lower bound we obtain lies within reach of the future generation 
of neutrinoless double beta decay experiments. 

The panels in Fig. [2] show the CP- violating phase S and the corresponding CP- violating 
invariant in neutrino oscillations: 

J = S12S23S13C12C23C13 sin (5 , (18) 



versus sin 2 #13. Note that these hold both for normal and inverted hierarchy spectra. 
In the middle panel one sees that cos (5 < since, at first order in sin 2 #13, m\jra 2 = 



1 + 



cos g 2 3 



cos 0i 9 sin 9-\ •> sin 8dh 



sin # 13 cos S, and the ratio of masses should satisfy: mi/m 2 < 1. 




FIG. 2: CP violating phase S and CP-invariant J in terms of the reactor mixing parameter. The 
3 o"-excluded range for sin 2 Oij is given for comparison. 



Moreover, for large 6* 13 values, where CP violation is likely to be probed in neutrino 
oscillations, one can see that our model predicts maximal violation of CP. Quantitatively, 
from the right panel one sees that the 3<r bound on #13: sin 2 #13 < 0.053 implies an upper 
bound: \J\ < 0.06 on the CP-invariant. 

In addition, the two-zero texture structure of our neutrino mass matrix may have other 
implications, for example for the expected pattern of lepton flavor violating decays. In 
fact, thanks to the strong renormalization effects due to the presence of the triplet states, 
the latter are quite sizeable in sypersymmetric models 
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V. MAJORON DARK MATTER 



In models where neutrinos acquire mass through spontaneous breaking of an ungauged 



lepton number 22|, |24j one expects that, due to non-perturbative effects, the Nambu- 



Goldstone boson (Majoron) may pick up a mass that we assume to lie in the kilovolt 



range 



40(. This implies that the Majorons will decay, mainly in neutrinos. As the coupling 



gj vv is proportional to — 24|, the corresponding mean lifetime can be extremely long, 
even longer than the age of the Universe. As a result the Majoron can, in principle, 
account for the observed cosmological dark matter (DM). 



This possibility was explored in Refs. 4l|, |42| in a general context. Here, we just 



summarize the results. It was found that the relic Majorons can account for the observed 
cosmological dark matter abundance provided 



T Jvv < 1.3 x 10~ 19 s- 1 , 0.12 keV <(3mj< 0.17 keV, 



(19) 



where Yj vv is the decay width of J — > vv and mj is the Majoron mass. The parameter (3 
encodes our ignorance about the number density of Majorons, being normalized to (3 = 1 



if the Majoron was in thermal equilibrium in the early Universe decoupling sufficiently 



early, when all other degrees of freedom of the standard model were excited 42| . In the 
following we will follow their choice and will take 

10" 5 < p < 1, (20) 

and calculate both the width into neutrinos as well as the subleading one-loop induced 
decay into photons. 

A. Decay into neutrinos 

We now proceed with the computation of the Majoron decay width into neutrinos, 
which will be useful to obtain the allowed parameter space for which the Majoron can be 
a viable DM candidate. In order to calculate the decay amplitude we remind that the 
coupling gj ViV is defined through 

C = —^gjuiVjJviVj + h.c. (21) 



For the evaluation of gj ViV , we follow the steps developed in Ref. [24(. First we 
notice that with scalar potential defined in section IHH the Majoron, in the basis 
[lm(0°), Im(A°), Im(cr )] , is given by 



J = Nj 
and 



2u A^ 2u A^ 2u \—pv U A^ 2 , U a (4u 2 A + V 2 ) 



(22) 



Nj = [4v 2 u A + v%l + ul{Aul + v 2 f] ~ 1/2 ~ -J- , (23) 

where the last equality follows from the assumed hierarchy ua «C v ^C u a implied by the 
vev-seesaw relation. Using this, one can obtain 

gju^ - — /= — , l^4j 

V 2 u a 
leading to the decay width 

ij ^"327r 2ul ■ [Zb) 

It is worth mentioning that the sum Yli( m i) 2 ^ s * n our framework constrained by the 
special form of the effective neutrino mass matrix shown in Eq. (TT6|) . In particular, there 
is a lower bound on the mass of the lightest neutrino: m > 0.03 eV. 



B. Decay into photons 



The Majoron also couples with photons (at the quantum level) and therefore the radia- 
tive decay J — > 77 is expected to occur with a photon energy E^ ~ mj/2. Consequently, 
this decay exhibits a mono-energetic emission line which could be detected in a variety of 
X-ray observatories, see for example the discussion given in Refs. 3l|, |42j. 

The effective Majoron-photon interaction can be written as 
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FIG. 3: Top: One loop diagrams for the decay J — > 77. Bottom: Effective J77 vertex. 



gjyryE 



± (mu- 1 a 



Pi 



(26) 



resulting from the one-loop diagrams shown in Fig. [3] (top diagrams). The effective cou- 
pling (yfj 77 (bottom graph in Fig. [3]) is 



9j-ry 



N f a 2 g Jff Q 2 f X f 



mm, 



(27) 



with Xf = —2m 2 rCo{0,0,m 2 j,m 2 r,m%,m 2 r) ~ 1 + m 2 /(12m 2 f) where Co is the invariant 
Passarino-Veltman loop function 43[. The last approximation is valid for mj <$; rrif. 
T3 , Qf and Nf denote the weak isospin, the electric charge and the colour factor of 
the corresponding charged fermionf , respectively. The coupling of the Majoron to the 
charged fermions gjff is given by 42J 



9Jff 



2u\ 

v 2 u a 



m f { 



-2T() 



(28) 



We then get for the decay width, 
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(29) 



where the cancellation of the anomalous contribution has been taken into account. 



Numerical results 



In this section we discuss some numerical results regarding the implementation of the 
decaying Majoron dark matter hypothesis in our scenario. In Ref. 42| it was shown that 
the experimental limit in the Majoron decay rate into photons is of the order of 1CT 30 s _1 . 
It was also shown that, in a generic seesaw model, a sizeable triplet vev plays a crucial 
role in bringing the decay rate close to this experimental bound. Here we have computed 
the width of the Majoron into neutrinos and photons in our extended seesaw model which 
incorporates the A4 flavor symmetry, generalizing the models of Ref. 29]. The results are 
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FIG. 4: Left panel: Tj w as function of the Majoron mass respecting Eq. (jT9j) for ua =1 eV 
(turquoise), 100 eV (dark green), lOkeV (magenta), lMeV (grey), lOMeV (dark blue) and 100 
MeV (black). Right panel: Tj 77 as function of the Majoron mass for the same values of the 
triplet vev as in the left panel. The upper orange shaded region is the excluded region from 
X-ray observations taken from Ref. 311 ] . 



shown in Fig. HI These take into account the current neutrino oscillation data, discussed 
in section IIV1 We chose five values for the triplet vev, u& =1 eV (turquoise), 100 eV 



(dark green), 10 keV (magenta), 1 MeV (grey) and 10 MeV (dark blue) and 100 MeV 
(black). For the right panel we consider only points that satisfy the WMAP constraint 
({T9l indicated by the red horizontal band on the top of the left plot. 

In order to be able to probe our decaying Majoron dark matter scenario through the 
mono-energetic emission line one must be close to the present experimental limits on 



the photon decay channel, discussed in Ref. [42| and references therein. As mentioned, 
this requires the triplet vev to be sizeable, as shown on the right panel of Fig. H]for the 
same choices of u&. In principle there is an additional lower bound on the Majoron mass 



coming from the Tremaine-Gunn argument 44j, which, for fermionic dark matter would 



be around 500 eV. Under certain assumptions this bound could be extended to bosons 



and is expected to be somewhat weaker 45j. The upper orange shaded region is the 



excluded region from X-ray observations given in Ref. 31] . One should point out that, in 
this model, because of the vev seesaw relation u&u a ~ v 2 one cannot arbitrarily take large 
values for «a to enhance Tj 77 because then the singlet vev gets correspondingly smaller 
values, hence reducing the lifetime of the Majoron to values in conflict with the WMAP 
constraint. This interplay between the CMB bounds and the detectability of the gamma 
line is illustrated in Fig. HJ where the dark-blue points corresponding to ua = 10 MeV 
illustrate the experimental sensitivity to our signal. 

VI. CONCLUSIONS 

We have studied the possibility that the seesaw model with spontaneously broken 
ungauged lepton number may simultaneously account for the observed neutrino masses 
and mixing as well as the dark matter of the Universe. We have presented a two-texture 
structure for the neutrino mass which arises in a specific seesaw scheme implementing an 
A4 flavor symmetry. A predictive pattern of neutrino masses emerges from the interplay 
of type-I and type-II seesaw contributions, with a lower bound on the neutrinoless double 
beta decay rate, which correlates with the deviation from maximality of the atmospheric 
mixing angle #23? as well as nearly maximal CP violation, correlated with the reactor 
angle 13 . 

On the other hand, assuming that associated Majoron picks up a mass due to explicit 
lepton number violating effects that may arise, say, from quantum gravity, we showed 
how it can constitute a viable candidate for decaying dark matter, consistent with cosmic 
microwave background lifetime constraints that follow from current WMAP observations. 
We have also shown how the Higgs boson triplet, whose existence is required by the 
consistency of the model, plays a key role in providing a test of the decaying Majoron 
dark matter hypothesis, implying the existence of a mono-energetic emission line which 
arises from the sub-leading one-loop-induced decay of the Majoron into photons. We 



also discussed the possibility of pro bing its existence in future X-ray observations such as 



expected in NASA's Xenia mission [46] . The presence of the type-II seesaw Higgs triplet 
would also have other particle physics implications, such as lepton flavor violating decay 
rate enhancements due to the strong renormalization effects of the triplet, quite sizeable 
in a supersymmetric model. 
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Appendix A: Basic A4 results 

The group A 4 consists of the even permutations of four elements and has three one- 



dimensional representations and one three-dimensional, see, e.g. J47J. Using the usual 
notation for transpositions and cyclic permutations (for instance, (123)4 applied to abed 
gives bead), the one-dimensional representations are shown in Table HT1 where ui = e 2nl ' 3 

TABLE II: Unidimensional representations for A4. 
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is the cubic root of unity, and the equivalence classes are defined as 



Class 1 
Class 2 
Class 3 
Class 4 



(12)(34), (13)(24), (14)(23) 
1(234), 2(143), 3(142), 4(132) 
1(243), 2(134), 3(124), 4(123). 



It follows immediately that 



V x 1' = 1", 1' x 1" = 1, 1" x 1" = 1'. 

As for the decomposition for the tensorial product of two triplets in A^ one has: 

3©3= l©l'©l"©3 s ©3 a , (Al) 
where the triplet and singlet representations are 

(u © v)\ = U 1 V 1 + U 2 V 2 + U 3 V 3 (A2) 

(u®v)v = U\V i + U 2 U 2 V2 + COU 3 V 3 (A3) 

(u®v)ii> = UiVx + luu 2 V2 + co 2 u 3 v 3 (A4) 

(«®v) 8s = (^2^3 + W3M2, M3W1 +^1^3, ^1^2 + U 2 Vi) (A5) 

(«<2)f) 3a = («2^3 - v 3 u 2 ,u 3 v 1 -v 1 u 3 ,u 1 v 2 - U 2 Vi). (A6) 
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